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Abstract 

We summarize some characteristic features of the frustrated magnetic interactions in 
spin-orbital models adequate for cubic transition metal oxides with orbital degeneracy. A 
generic tendency towards dimerization, found already in the degenerate Hubbard model, 
is confirmed for tig but not for e g systems. In the t2 g case the quantum orbital fluctu- 
ations are more pronounced and contribute to a stronger competition between different 
magnetic and orbital states. Therefore the orbital liquid states exist in some undoped 
t2 g systems, while in the manganites such states can be triggered only by doping. 
Journal reference: A. M. Oles, Phys. Stat. Sol. (b) 236, 281 (2003). 

The physical properties of transition metal oxides are dominated by large on-site Coulomb 
interactions oc U which suppress charge fluctuations. Therefore, such systems are either Mott 
or charge-transfer insulators, and the metallic behavior might occur only as a consequence 
of doping. Here we will discuss first the undoped systems with localized d electrons which 
interact by effective superexchange (SE) interactions. An interesting situation occurs when 
d electrons occupy partly degenerate orbital states, and one has to consider orbital degrees 
of freedom in the SE at equal footing with electron spins [1]. Competition between different 
states is then possible, holes may couple to orbital excitations [2], and the quantum effects 
are enhanced already in undoped systems [3] . The first models of SE in such situations were 
proposed almost three decades ago [4], either by considering the degenerate Hubbard model 
[5, 6], or for realistic situations encountered in cuprates (KCUF3 and K2CUF4) and in V2O3 
[7]. Then it was realized that the SE which is usually antiferromagnetic (AF) might become 
ferromagnetic (FM) when Hund's exchange interaction Jh is finite, but only in recent years 
the phenomena which originate from the orbital physics are investigated in a more systematic 
way. 

The SE which involves the orbital degrees of freedom is described by the so-called spin- 
orbital models [8], and is typically highly frustrated on a cubic lattice where it might even 
lead to the collapse of magnetic long-range order by strong spin or orbital fluctuations [3]. 
However, in real e g systems such quantum phenomena are usually quenched by finite Jh 
which induces a structural phase transition and thus helps to stabilize a particular ordering 
of occupied orbitals which supports A-type AF order, as observed when degenerate orbitals 
are filled either by one hole (KCUF3) [9], or by one electron (LaMnOs) [10]. The coupling to 
the lattice due to the Jahn- Teller (JT) effect also helps to stabilize the orbital ordering, and 
quantitative models of the structural transition have to include both these effects [10]. 

The essential feature of the SE described by spin-orbital models is the frustration of mag- 
netic interactions: the FM terms occur next to the AF ones, and it depends on the physical 



1 ) E-mail: A.M.Oles@fkf.mpg.de 



phys. stat. sol. 



- 2 - 




i 

U+2J H 




'2 


U 
U-J H 


: 1 i*L e 


2 liL E 


U-3J H 








d 2 


d 3 



Figure 1: Excitation spectra in cubic transition metal oxides for: (a) e g systems: Cu 3 "*" 
(d 8 ) and Mn 2 + (d 5 ) ions; (b) t 2g systems: Ti 2 + (d 2 ) and V 2 + (d s ) ions [8]. 



parameters which interactions finally win and stabilize a given type of magnetic order. The 
simplest spin-orbital model which illustrates this physics may be derived for the density of 
one electron per site (n = 1) in a doubly degenerate Hubbard model. We assume that the 
hopping is isotropic and diagonal between two orbitals a and /3 at sites i and j, and include 
only the on-site interaction elements: Coulomb U and Hund's exchange Jh, 

H = -t ^a-^a^CT + ?7^n iMT n if4 + (U - ^Jff)^n iQ n 4(3 

- 2J H S la Sn3 + J H ^(aL T aL| a i/3| ffl i/3T + a \^ a \i3i a iai a ia]), (1) 



with the spin, {St, S ifJ , Sf M } = {al^a^, a^^^ , (n <MT - n iMi )/2}, and density operators, 
nip = n i^] + n iMii a ^ orbital fi — a, (3 defined in the usual way. The interactions are 
rotationally invariant in the orbital space [11]. If U S> t, the electrons localize and the low- 
energy physics is described by the SE interactions which follow from the virtual djdj ^ d ( 'd 2 
processes on the bonds (ij). They lead either a high-spin 3 A2 state, or to a double occupancy 
in either orbital which has to be subsequently projected onto two low-spin 1 E and 1 A\ 
states. The excitation spectum is equidistant, as for d ions in the cuprates [Fig. 1(a)], 
with the excitation energies derived from Eq. (1): £(^2) = U — 3Jh, s(^E) = U — Jh, and 
e(^4i) — U + Jh [9]. The SE Hamiltonian derived from Eq. (2) takes the form, 



Hi 



= Jr 1 ^(s i -S j + l)(f i .f j -\) + 1 -J^(s i .S j -\) 

(ij) (ij) 

x \r 2 {l+T+T-+T-T+) + r 3 (i + 27? T/) + (r 2 - r 3 )(T?T? +TrTf) 



,(2) 



where J = 4t 4 /U is the energy unit for the SE interaction, and the coefficients n = 1/(1 — 377), 
T2 = 1/(1—77), r 3 = 1/(1+77) follow from the above charge excitations, where 77 = Jh/U. Simi- 



i0ai a 'i{3 a a ia<7, i n i' 



lar to spin, the pseudospin operators are: {T^ + , T~ , If } = {^<t a iarT a i0ai ^2 
7ij / g)/2}. It is important to use the accurate form of the electron-electron interactions [5, 9], 
and for this reason some early work led to inaccurate expressions [4, 6]. Note that spin 
interactions have SU(2) symmetry, while the orbital interactions are anisotropic. The first 
term is simple and follows from the excitations of spin triplet and interorbital singlet state. 
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Figure 2: (a) Mean-field phase diagram in (Jh,Eq) plane of the isotropic spin-orbital 
model (2) with four different magnetic phases: DIM, FM, AF1 and AF2 (with cither a or 
P orbitals occupied); (b) exchange constants J a (, and J c for increasing Jh at Eo = as 
obtained in the DIM phase (Jfm and Jaf) and in the FM phase (J ab = J c ). 



The low-spin terms oc 7*2(3) are more involved and include not only orbital-flip processes, but 
also pair hopping terms oc (T^T^ +T~TJ~). This demonstrates that the anisotropy in the 
orbital sector is a feature which follows from the multiplet spectra of transition metal ions 
[12], where the orbital triplet state never occurs at Jh > 0. 

The model (2) simplifies in the limit of Jh — > 0, and represents a superposition of excita- 
tions which involve either spin triplet and orbital singlet, or spin singlet and orbital triplet, 



ff,=j£ [(&.$,■+!) (£.3 -i) -i) {f i .f j + 

(ij) 



(3) 



which is just a different way of writing the SU(4) symmetric spin-orbital model [13]. In 
this case the spin and orbital correlations obey full SU(4) symmetry, and the correlations 
functions: (Si ■ Sj), (f; ■ fj), §((S; • Sj)(fi ■ fj)), are all identical [14]. This condition is 
violated when the mean-field approximation (MFA) is used and the spin and orbital variables 
are decoupled, so the results of the MFA might be unreliable. 

Although some qualitative arguments were given, the classical phase diagram of the spin- 
orbital model (2) was not investigated before. We include the orbital splitting at every site, 
~ Eo(n ia —mp)/2, and compare the energies of four different three-dimensional (3D) phases: 
(i) AF long-range order (LRO) with either a or /3 orbital occupied at every site, (ii) FM phase 
with alternating a//3 orbitals on two sublattices, and (Hi) a dimer phase (DIM) characterized 
by orbital valence bond (OVB) states, with orbital singlets at every second bond along c axis 
(or any other, as the present problem is isotropic). When the orbital singlet is formed on 
a single bond, the energy gain due to the first term in Eq. (2) is maximized and the FM 
interaction follows. This leads at small Jh to a DIM state in the model for vanadates [15], 
as we will discuss below. On the contrary, the orbitals are uncorrelated at all other bonds, 
the AF terms win as long as the Hund's interaction is weak (Fig. 2). The AF states: AF1 
and AF2 are stabilized by the orbital splitting Eo which has to counterbalance the energy 
gains on the FM bonds. Of course, it is hard to imagine that the DIM state with ordered 
orbital singlets might be realized as such, but the alternation of FM/AF bonds is plausible, 
so the present phase diagram should rather be viewed as demonstrating a generic competition 
between different signs of the SE interactions. It shows that one may indeed expect enhanced 
quantum fluctuations close to the orbital degeneracy when Jh is small [3] . 

The simplest realistic spin-orbital model describes d 9 ions interacting on a cubic lattice, 
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as in KCUF3. The interactions are the same as in Eq. (1), but the hopping term is now 
nondiagonal and allows for orbital excitations [2] . In the limit of U 3> t the charge excitations 
dfdj ^ dfd)° lead again to the same excited states as above, with their energies [Fig. 1(a)] 
reproducing the exact spectrum of d ions [12]. We define the SE J = 4t 2 /U by the largest 
hopping element t between two \z) = \3z 2 — r 2 ) orbitals along the c axis, and one finds, 

^ 9 ) = i J E E [(&-3 + i)W)+*£V)]> (4) 

7 (ij>ll7 

where 7 = a,b, c, and Si are spin S = 1/2 operators. The operator expressions: 

J<7V) = (2 + » 7 r2-V3)P < %-i;(3ri-r 2 )7'« ) , (5) 

*&V) = -[l + »/(3ri+r 2 )/2]P« - [1 + „(r 2 - r 3 )/2]P« , (6) 

describe spin and orbital SE, and the coefficients n are defined as in Eq. (2). The operators: 

V% } = (l/2 + r7)(l/2-r7) + (l/2-r7)(l/2 + r7), (7) 

7>« = 2(1/2 -r7)(l/2-r/), (8) 

project on the orbital states, being either parallel to the bond (ij) direction on one site 
(Pit = 1/2 — r t 7 ) and perpendicular on the other (Pj£ = 1/2 + tJ) one, or parallel on both 
sites. They are represented by the orbital operators r 7 for the three cubic axes: 



i(b) 



(-a*±V3<rf)/4, T f = o* /2, (9) 



where the cr's are Pauli matrices acting on: |a;) = ^ J ^ , \z) = ^ J ^ 1 which transform as 

|x) tx x 2 - and |z) oc (3z 2 - r 2 )/\/3. 

In LaMnOa the SE is more involved and couples total spins S — 2 at the Mn 3+ ions. 
It originates from the charge excitations, dfdlj ^ dfd'j [10]. The e s part, following from 
dfdj ^ df (^2g)<^i (^2 9 e g) processes, involves again FM terms due to the high-spin r 'Ai state, 
and three AF terms due to the low-spin states: 4 Ai, 4 E, and 4 A2 [Fig. 1(a)], and has 
analogous orbital dependence as in the cuprate case. In contrast, the t2 g part follows only 
from low-spin excitations d 4 d 4 ^ df (t2 g )d^(t 4 g e g ) and is therefore AF and almost orbital 
independent. Both terms are given explicitly in Ref. [10]. 

Both the cuprate model (4) and the e g term in the manganite model describe strongly 
frustrated SE interactions, which take a universal form in the limit of Jh — * 0, 



"8* -i's; 0(5-*) (5 



(10) 



1 <ij> IIt 



Several classical phases have the same energy of —3 J/4 per site [3]: the G-AF phases with 
arbitrary occupation of orbitals, and A-AF phases with ((1/2— t?)(1/2— tJ)) = 0, as obtained 
for staggered planar orbitals, e.g. for x 2 — y 2 /y 2 — z 2 orbitals staggered in (a,b) planes. We 
emphasize that the model (10) is qualitatively different from the idealized SU(4) symmetric 
case (3) due to the directionality of e g orbitals. In fact, the e g orbitals order easier, may 
couple to the lattice and thus appear to be more classical than the isotropic case described 
by Eq. (2). Their ordering supports magnetic phases with coexisting FM [in (a, b) planes] 
and AF (along c axis) interactions. 

The classical phase diagram of the cuprate model (4) is shown in Fig. 3(a). Quantum 
corrections to this phase diagram were discussed in Ref. [3]. They suggest that a spin liquid, 
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Figure 3: Mean-field phase diagrams for: (a) one-hole model in (Jh,E z ) plane, and 
(b) one-electron d 1 model in (Jh,E c ) plane. In the cuprates (a) the magnetic interactions 
are predominantly AF and the order changes from a 2D AFrr through A-AF to a 3D AFz 
phase, while in the titanatcs (b) the DIM, FM, and AFc phases are stable. 



supported by particular OVB correlations, might be realized near the degeneracy of classical 
phases. Finite Jh stabilizes the A-AF phase, with staggered two-sublattice orbital order, 
\incr) = cos6i\iza) ± sm6i\ixa) , where ± refers to i € A(B) sublattice. The AF interactions 
decrease with increasing Jh/U and dominate at realistic Jh/U ~ 0.12 [18]. However, the 
FM interactions within the (a, b) planes predicted by the model [8] are considerably stronger 
than measured [16], showing that the quantitative understanding requires also Goodenough 
processes which would provide additional AF interactions. In contrast, the FM interactions 
are stronger than AF ones in the A-AF phase realized in LaMn03 [17], and are much better 
reproduced by the SE terms derived in Ref. [10] with rj = 0.117 [18]. 

The transition metal oxides with partly filled t2 g orbitals exhibit different and even more 
interesting phenomena. In this case the JT coupling is much weaker, and (unlike for the e 9 
orbitals) the orbital quantum number is conserved in the hopping processes. This leads to 
qualitatively different physics realized in t2 g systems, somewhat more similar to the isotropic 
case, Eq. (2). Each t2 g orbital is orthogonal to one of the cubic axes, so we label them as a, 
b , and c (for instance, xy orbitals are labelled as c) . The models for titanates and vanadates 
follow from the d"d" =^ d™~ 1 d™ +1 processes [19, 20], and may be written in a general form: 



7 (ij)\h 



(11) 

with the exchange constants J^\d n ) between S = 1/2 spins for titanates (n = 1) and 



5=1 spins for vanadates (n — 2), and purely orbital interactions K>J (d n ). In titanates 
these interactions depend on the Hund's rule splittings of d 2 ions [12] [Fig. 1(b)] via the 
coefficients: n = 1 — 3r/, — 1 — rj, rs = 1 + 2i], and were given in Refs. [8, 21]. They 
are faithfully reproduced with a model Hamiltonian containing U and Jh for t2 g orbitals. 
A priori, the magnetic interactions are anisotropic, and may be either AF or FM, depending 
on the orbital correlations. 

In the limit of Jh/U = the Hamiltonian (11) takes the form, 



(12) 



7 (»J>||7 



and shows again a strong frustration of SE interactions [19]. Although it resembles formally 
the SU(4)-symmetric spin-orbital models [13] even more than Eq. (10), the pseudospin oper- 
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Figure 4: (a) Mean-field phase diagram in (Jh,E c ) plane as obtained for the vanadate 
model of Ref. [20]. (b) Exchange constants for increasing Jh at E c = for the DIM (J c i, 
Jc2, J a b) an d f° r C-AF (J c and J a f,) phase [27]. 

ators r» = {rf , t? , t*} have here a different meaning and refer to a pair of orbital flavors for 
each cubic direction 7, given by two active t2 g orbitals which contribute to the SE [19, 20]. 
Thus, the model is again different from the idealized SU(4) symmetry Eq. (3). 

At finite rj the magnetic interactions are a priori anisotropic, and may be either AF or 
FM, depending on the orbital correlations. In order to get some qualitative insight into the 
competition between these terms, one may include an anisotropy term ~ E c [rii c — (n ia + 
n,b)/2], and evaluate the energy of a DIM phase with fluctuating a and b orbitals along c 
axis (riia + nn, — 1), a FM phase with equally and randomly occupied orbitals (n; 7 = 1/3), 
and a two-dimensional (2D) AF phase with only c orbitals occupied (m c = 1). Unlike for d 9 
case, the FM phase is stable in a broad regime of parameters [Fig. 3(b)], and the DIM phase 
is stabilized by E c ~ 0.2J. Of course, this analysis is oversimplified and large corrections 
due to quantum effects are expected. Indeed, a FM isotropic phase is realized in YTi03, 
but a closer inspection shows that the orbitals do order, but this ordering does not break 
the cubic symmetry [22]. A completely different state is realized in LaTiC>3, however, with 
isotropic AF interactions [23]; such interactions are explained by quantum resonance realized 
simultaneously in spin and orbital sector [19]. This shows that a particular type of magnetic 
ordering may be triggered by quantum fluctuations in the orbital liquid. 

The magnetic ordering realized in vanadates is different: C-type of AF order is observed 
both in LaVC>3 [24] and in YVO3 at intermediate temperatures 77 < T < 116 K, and G-type 
AF order is stable in YVO3 for T < 77 [25]. As in V2O3 [26], the SE interactions between 
5=1 spins follow from the dfd| =^ d\d 3 j processes, leading to the effective spin-orbital model 
given by Eq. (11) with n = 2. When the electrons condense in c orbitals (n» c = 1) due to the 
orbital splitting caused by the JT effect, the second electron occupies either a or 6 orbital at 
every site (n ia +riib = 1), allowing for a resonance on the bonds (ij) along c axis. In this case 
the pseudospin operators are: rf = ajbj, rf = bjcij, t* = \{n ia —n ib ), and n' c ' = n ia + n ib , 
where {a\,b\} are Schwinger bosons for a and 6 orbitals at site i. 

The vanadate model has again an interesting classical phase diagram which unifies certain 
features we have already seen for degenerate isotropic orbitals, and in the t\ g model [Fig. 
4(a)]. At r) — and E c — one finds again the frustrated SE (12) between S = 1 spins. 
While the orbital liquid cannot stabilize in this case, orbital singlets may form along the c 
direction when c orbitals have condensed (n; c = 1) and the a and 6 orbitals fluctuate. This 
favors the OVB state, with strong FM interactions alternating with weak AF ones along the 
one- dimensional (ID) chains [15]. At large Jh this state is unstable, however, and the orbital 



phys. stat. sol. 



-7- 




X 



Figure 5: (a) Spin-wave stiffness D as a function of hole doping x obtained from the DE 
(dashed line), and by including also the SE terms (solid line) for: A-AF (AF), FI, and 
FM metallic phase [30]. Experimental points for: Lai-^Sr^MnOa [31] (diamonds) and 
Lao.7Pbo.3Mn03 [32] (filled circle), (b) Magnon dispersion ui^ for x = 0.30 (solid line) 
[30], and the data points for Lao.7Pbo.3Mn03 (circles and dashed line) [32]. 

fluctuations support FM interactions along c axis and stabilize the C-AF phase [20]. By 
considering the energy in the MFA one finds a phase transition from the DIM to C-AF phase 
at rj c ~ 0.09 [27], as long as rii C — 1. At large uniform orbital splitting E c > 0, the charge 
gets redistributed to n; a = mb = 1, and an anisotropic G-AF state with the strong AF bonds 
along c axis, and weaker ones within (a, b) planes, follows. To our knowledge, such a state 
has not been observed so far. The G-AF phase found in YVO3 at T < Tni is characterized 
by large JT distortions, and thus can be explained by a staggered field which favors G-type 
orbital ordering [20], in agreement with recent experiments [25]. 

The SE interactions in transition metal oxides depend on the multiplet splittings oc Jh- 
For example, if n ic — 1 and n ia + mb = 1 in cubic vanadates, the exchange constants within 
the (a, b) planes (J a t) and along the c axis (J c ) in the G-AF phase are given by [20]: 

Jab = [l-rj(R+r) + {l + 2rjR-rjr){n la n ja ) {b) ]/4:, (13) 
J c = [(1+2^X7^ + 1/4)^- W (t?t* + 1/4)^-7^/2, (14) 

where R — 1/(1 — 377) and r — 1/(1 + 2rf). Similar expressions can be derived in the 
DIM phase at r\ < r/ c . Assuming orbital singlets in the DIM phase, and strong ID orbital 
fluctuations, described by a pseudospin ID Heisenberg model in the G-AF phase, one finds 
that the FM/AF exchange constants coexist and increase/decrease with increasing Jh [Fig. 
4(b)]. It is interesting to observe that the values of J a b and \J C \ arc similar for a realistic 
value of r\ ~ 0.116 [18], as the orbital fluctuations enhance the FM interactions oc J c . 

Let us come back to the question why the orbital liquid state cannot stabilize in LaMn03 . 
In this case the orbitals do order, and the orbital interactions are so strong that their ordering 
would occur well above Tn even in the absence of the JT interaction [10]. However, the 
splitting between the high-spin 6 Ai state and low-spin states is 5 Jh (Fig- la), which explains 
the proximity to the FM ordering. The manganites at x < 0.15 are insulating [17], and are 
orbital ordered, with either ^4-AF or FM insulating (FI) phase due to polaronic effects [Fig. 
5(a)]. Therefore, a single hole does not propagate freely but scatters on orbital [28] and spin 
[29] excitations. We discussed elsewhere that doping x > 0.15 stabilizes the FM metallic 
state due to the double exchange (DE) for strongly correlated e g orbitals [30]. This FM 
metallic state is nothing else than the realization of the orbital liquid in an e g system. By 
considering the DE and SE together one arrives at a quantitative explanation of: (i) the spin- 
wave stiffness D increasing with x [31], which just reflects the gradual release of the kinetic 



phys. stat. sol. 



- 8- 



energy by hole doping [Fig. 5(b)]; (ii) the isotropic spin waves observed around x = 0.3 in a 
system being so susceptible towards the orbital ordering. The spin- wave stiffness D e g — 7.45 
meV obtained at x = 0.3 without any fitting parameters agrees well with D cxp = 8.79 meV 
measured in Lao.7Pbo.3Mn03 [32] and explains the observed dispersion tu^ [Fig. 5(b)]. A 
transition from the FM to A-AF phase observed in bilayer manganites La2-2aiSri + 2a;Mn207 
[33] can also be explained within the same approach. 

In summary, the transition metal oxides with orbital degrees of freedom show a very 
fascinating behavior, with various types of magnetic and orbital order. While e g orbitals 
usually order and explain A-AF phases, further stabilized by the JT effect, the tig orbitals 
have a generic tendency towards disorder, which leads to the isotropic orbital liquid in the G- 
AF phase in LaTiOa, and to a ID anisotropic orbital liquid in the C-AF phase in LaV03 and 
YVO3. So strong orbital fluctuations in e 9 systems and the orbital liquid state are triggered 
only by large doping in the manganites. Very interesting quantum effects might also soon be 
discovered in the orbital liquid states in doped titanates and vanadates. 

Acknowledgments. It is a pleasure to thank Lou-Fe' Feiner, Peter Horsch, Giniyat Khal- 
iullin, and Jan Zaanen for a very friendly collaboration on this subject and for numerous 
stimulating discussions. I thank also J. Bala, A. Fujimori, B. Keimer, G.A. Sawatzky, and 
Y. Tokura for valuable discussions. This work was supported by the Committee of Scientific 
Research (KBN) of Poland, Project No. 2 P03B 055 20. 

References 

[1] Y. TOKURA and N. NAGAOSA, Science 288, 462 (2000). 
[2] J. Zaanen and A. M. Oles, Phys. Rev. B 48, 7197 (1993). 
[3] L. F. Feiner, A. M. Oles, and J. Zaanen, Phys. Rev. Lett. 78, 2799 (1997). 
[4] K. I. KUGEL and D. I. Khomskii, Sov. Phys. JETP 37, 725 (1973); Sov. Phys. Usp. 25, 231 
(1982). 

[5] M. Cyrot and C. Lyon-Caen, J. Phys. (Paris) 36, 253 (1975). 
[6] S. INAGAKI, J. Phys. Soc. Jpn. 39, 596 (1975). 

[7] C. Castellani, C. R. Natoli, and J. Ranninger, Phys. Rev. B 18, 4945, 4967 and 5001 (1978). 
[8] A. M. Oles, Acta Phys. Polon. B 32, 3303 (2001). 

[9] A. M. Oles, L. F. Feiner, and J. Zaanen, Phys. Rev. B 61, 6257 (2000). 
[10] L. F. Feiner and A. M. Oles, Phys. Rev. B 59, 3295 (1999). 
[11] A. M. Oles, Phys. Rev. B 28, 327 (1983). 

[12] J. S. Griffith, The Theory of Transition Metal Ions (Cambridge University Press, Cambridge, 
1971). 

[13] Y. Q. Li, M. Ma, D. N. Shi, and F. C. Zhang, Phys. Rev. Lett. 81, 3527 (1998). 

[14] B. Frischmuth, F. Mila, and M. Troyer, Phys. Rev. Lett. 82, 835 (1999). 

[15] S. Q. Shen, X. C. Xie, and F. C. Zhang, Phys. Rev. Lett. 88, 027201 (2002). 

[16] B. Lake, D. A. Tennant, and S. E. Nagler, Phys. Rev. Lett. 85, 832 (2000). 

[17] F. Moussa, M. Hennion, C. Biotteau, J. Rodri'guez-Carvajal, L. Pinsard, and A. 

Revcolevschi, Phys. Rev. B 60, 12299 (1999). 
[18] T. Mizokawa and A. Fujimori, Phys. Rev. B 54, 5368 (1996). 
[19] G. Khaliullin and S. Maekawa, Phys. Rev. Lett. 85, 3950 (2000). 
[20] G. Khaliullin, P. Horsch, and A. M. Oles, Phys. Rev. Lett. 86, 3879 (2001). 
[21] G. Khaliullin, Phys. Rev. B 64, 212405 (2001). 

[22] G. Khaliullin and S. Okamoto, Phys. Rev. Lett. 89, 167201 (2002). 

[23] B. Keimer, D. Casa, A. Ivanov, J. W. Lynn, M.v. Zimmermann, J. P. Hill, D. Gibbs, Y. 

TAGUCHI, and Y. Tokura, Phys. Rev. Lett. 85, 3946 (2000). 
[24] S. Miyasaka, T. Okuda, and Y. Tokura, Phys. Rev. Lett. 85, 5388 (2000). 
[25] Y. Ren, T. T. M. Palstra, D. I. Khomskii, A. A. Nugroho, A. A. Menovsky, and G. A. 

Sawatzky, Phys. Rev. B 62, 6577 (2000). 
[26] R. Shiina, F. Mila, F.-C. Zhang, and T. M. Rice, Phys. Rev. B 63, 144422 (2001); S. Di 

Matteo, N. B. Perkins, and C. R. Natoli, Phys. Rev. B 65, 054413 (2002). 



phys. stat. sol. 



- 9- 



[27] A. M. Oles, P. Horsch, and G. Khaliullin, Acta Phys. Polon. B 34, 857 (2003). 

[28] J. van den Brink, P. Horsch, and A. M. Oles, Phys. Rev. Lett. 85, 5174 (2000). 

[29] J. Bala, G. A. Sawatzky, A. M. Oles, and A. Macridin, Phys. Rev. Lett. 87, 067204 (2001). 

[30] A. M. Oles and L. F. Feiner, Phys. Rev. B 65, 052414 (2002). 

[31] Y. Endoh and K. Hirota, J. Phys. Soc. Jpn. 66, 2264 (1997). 

[32] T. G. Perring, G. Aeppli, S. M. Hayden, S. A. Carter, J. P. Remeika, and S.-W. Cheong, 

Phys. Rev. Lett. 77, 711 (1996). 
[33] T. G. Perring, D. T. Adroja, G. Chamboussant, G. Aeppli, T. Kimura, and Y. Tokura, 

Phys. Rev. Lett. 87, 217201 (2001). 



